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1.0  INTRODUCTION 


This  report  considers  the  use  of  aggregation  in  modeling  and  solving 
network  based  logistics  planning  problems.  Aggregation  combines,  or 
pools,  certain  attributes  of  a  problem  to  facilitate  its  modeling  and 
solution.  Frequently,  after  the  model  is  solved  it  is  necessary  to 
disaggregate  the  pooled  attributes  to  their  original  level  of  detail. 

This  allows  determination  of  the  acceptability  of  the  model  output  in 
making  the  suggested  decisions.  If  the  aggregate  model  yields  an 
unacceptable  solution,  a  mechanism  is  required  by  which  the  solution  can 
be  modified  to  make  it  acceptable. 

Aggregation  imposes  a  natural  decision  hierarchy  with  higher  level 
aggregate  decisions  restricting  lower  level,  more  detailed,  decisions. 

For  large  scale  planning  situations  requiring  involvement  of  a  number  of 
participants  with  skills  in  different  areas,  this  can  be  very  important. 

A  structure  is  required  to  coordinate  the  efforts  of  the  different 
participants.  This  structure  should 

-  provide  initial  guidelines  for  aggregation, 

-  utilize  detailed  plans  for  revising  the  aggregate  plan, 

-  iterate  until  an  acceptable  overall  plan  is  obtained. 

A  mathematical  structure  for  the  aggregation  process  will  be  developed 
which  allows  a  large  problem  to  be  solved  as  a  sequence  of  smaller 
problems.  This  structure  can  also  serve  as  a  coordinating  instrument  for 
the  decision  making  hierarchy. 

Motivation  for  examining  aggregation  in  the  context  of  network  based 


planning  problems  is  summarized  as  follows.  Detailed  data  is  often  not 
available  at  the  beginning  of  the  planning  process,  but  aggregate  data 
can  be  estimated.  In  such  cases  solving  the  aggregate  problem  frequently 


Indicates  what  detailed  data  Is  required  before  a  solution  can  be 
implemented. 

Aggregating  data  tends  to  reduce  variation.  For  example,  monthly 
demand  is  likely  to  have  less  variation  than  dally  demand.  Consideration 
of  a  problem  at  a  detailed  level  often  requires  that  parameters  be 
modeled  as  random  variables,  while  an  aggregate  level  allows  the 
parameter  to  be  considered  as  essentially  deterministic .  This 
significantly  decreases  model  complexity. 

In  many  situations  a  detailed  planning  model  is  too  large  to  be 
solved  with  available  computer  resources  and  planning  time.  Further,  in 
many  planning  environments  the  objective  is  not  merely  to  optimize  a 
given  function  but  to  generate  a  plan  which  is  satisfactory  with  respect 
to  multiple  objectives,  some  of  which  are  not  quantitfiable .  All  or 
parts  of  the  model  may  be  solved  more  than  once,  with  the  user 
interacting  to  guide  the  process  to  an  acceptable  solution.  This  implies 
a  need  for  reorganizing  a  large  monolithic  model  into  a  set  of  smaller 
models  which  when  linked  together  represent  the  original  model. 
Aggregation  can  be  used  to  devise  such  a  decomposition  strategy. 

Aggregation  concepts  developed  here  are  motivated  by  the  need  for 
aggregation  In  the  SCOPE  system  developed  In  FDRC  Report  84-09.  One  of 
the  fundamental  modeling  components  of  the  SCOPE  system  Ip  the  assignment 
of  movement  requirements  to  transportation  channels  given  the  MR 
assignment  problem.  The  optimization  model  used  is  a  transportation 
model.  This  model  will  be  referred  to  throughout  this  report.  It  is 
discussed  in  detail  in  PDRC  Report  84-09.  The  generalized  version  of 
this  model  is  currently  being  coded  as  the  MRMATE  model  in  the  MODES 
system  under  developed  for  the  Joint  Deployment  Agency.  The  methodology 


developed  here  Is  applicable  to  both  the  pure  and  the  generalized  forms 
of  the  model  in  that  report.  To  simplify  the  presentation  here,  standard 
notation  will  be  used  rather  than  the  more  complex  notation  required  to 
describe  the  movement  requirement  assignment  model. 
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2.0  AGGREGATION 


Consider  the  example  transportation  problem  given  below  and  Its 
corresponding  network  representation  given  in  Figure  2.1. 
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In  the  model,  x^  is  the  flow  on  the  arc  from  source  node  1  to  sink 
node  j  and  c^  is  the  corresponding  cost.  Nodes  on  the  left-hand  side  of 
the  transportation  network  are  source  nodes  and  nodes  on  the  right-hand 
side  are  sink  nodes. 

This  report  concentrates  only  on  aggregation  of  nodes  on  the 
right-hand  side  of  the  transportation  network  (i.e.  sink  nodes  only). 

This  same  logic  can  be  extended  to  aggregation  of  both  sources  and  sinks 
by  performing  the  aggregation  in  two  steps. 

Consider  aggregation  of  nodes  4  and  5  in  Figure  2.1  to  form  an 
aggregate  node  1  in  Figure  2.2,  and  nodes  6  and  7  in  Figure  2.1  to  form 


aggregate  node  2  In  Figure  2.2.  Constraints  corresponding  to  the  network 
In  Figure  2.2  are  shown  in  (II)  below. 
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In  the  model,  y^  represents  the  flow  on  the  arc  from  source  node  i 
to  aggregate  sink  node  r  and  c£  is  the  corresponding  cost  per  unit  of 
flow. 

Examining  the  information  that  is  lost  in  the  aggregate  problem  in 
Figure  2.2  the  following  observations  can  be  made.  There  is  a  single 
cost  on  each  aggregate  arc  in  Figure  2.2.  This  implies  that  flow  on  the 
aggregate  arc  (1,1)  incurs  the  same  cost  as  it  would  on  either  of  the 
detailed  arcs  (1,4)  or  (1,5).  If  the  two  detailed  arc  costs  are  not  the 
same,  then  some  amount  of  error  is  introduced  by  combining  the  two  costs 
into  a  single  cost.  If  the  amount  of  flow  on  each  of  these  arcs  is  known 
in  advance,  the  aggregate  cost  could  be  weighted  so  that  it  would 
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accurately  reflect  the  detailed  cost.  While  the  detailed  flows  are  not 
known  in  advance,  there  is  the  potential  for  iteratively  improving  the 
weighted  aggregate  cost  after  solving  an  aggregate  problem  and  then 
disaggregating. 

A  simple  method  for  specifying  the  capacity  of  aggregate  nodes  is  to 
set  the  aggregate  capacity  equal  to  the  sun  of  the  corresponding  detailed 
node  capacities.  However,  such  a  capacity  may  not  be  meaningful  when  the 
nodes  model  different  entities.  For  example,  aggregating  sea  and  air 
channel  capabilities  in  the  MR  assignment  model  allows  materiel  to  be 
sent  by  either  air  or  sea.  If  some  materiel  can  only  go  by  sea  then  the 
aggregate  node  capacity  exceeds  the  practical  total  which  can  be  shipped. 

There  is  also  a  problem  in  specifying  the  capacity  of  aggregate 
arcs.  For  example,  consider  the  aggregate  arcs  (1,1)  and  (2,2)  in  Figure 
2.2.  Since  (1,1)  represents  both  arcs  (1,4)  and  (1,5)  in  Figure  2.1,  it 
can  have  a  maximum  flow  of  (a^  +  a^).  However  (2,1)  in  Figure  2.2 
represents  only  (2,5)  in  Figure  2.1.  Hence,  it  can  only  have  a  maximum 
flow  of  a^.  Procedures  must  be  developed  for  passing  this  information  to 
the  aggregate  problem  while  retaining  its  transportation  structure. 

Simply  putting  the  sun  of  the  associated  node  capacities  as  capacities  on 
the  aggregated  arcs  will  not  suffice. 

This  report  outlines  approaches  to  deal  with  each  the  above  issues. 

It  also  proposes  a  format  which  can  be  used  to  examine  alternative 
formulations  of  the  aggregate  problem  to  determine  an  optimum 
aggregation. 


/.  /v 
'  ./  ./ 
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2 . 1  General  Aggregation  Model 


This  section  describes  approaches  to  aggregation  methods  presented 
In  the  literature.  The  formulations  are  discussed  in  Zipkin  [5]  and 
Taylor  [A]. 

When  the  aggregate  problem  is  formulated,  the  following  parameters 
must  be  specified 

-  capacities  on  the  aggregate  sink  nodes 

-  costs  on  the  aggregate  arcs 

-  capacities  (upper  limits  on  flow)  on  the  aggregate  arcs. 

The  detailed  problem  is 
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Construction  of  the  aggregate  problem  is  essentially  a  two  step 
procedure.  First,  the  nodes  to  be  aggregated  (e.g.  Figure  2.3)  are 
combined  into  a  single  aggregate  node  (e.g..  Figure  2. A).  This  results 
in  a  set  of  parallel  arcs  as  indicated  in  Figure  2. A.  These  parallel 
arcs  are  then  aggregated  into  a  single  aggregate  arc. 

The  aggregate  node  capacity  is  derived  from  the  capacities  of  the 
individual  nodes  forming  the  cluster.  If  a  linear  relationship  is 
assumed  (as  a  simplifying  assumption),  then  the  aggregate  capacity  can 
expressed  as 


dip  Cj  bJ 


where  Pf  is  the  set  of  nodes  in  the  rth  aggregation  cluster,  is  the 
capacity  of  node  j  in  for  the  detailed  problem,  and  t^  is  a  user 
defined  multiplier  between  zero  and  one. 


The  motivation  for  different  choices  of  t^  are  as  follows.  If  there 
are  arcs  from  all  source  nodes  to  all  sink  nodes  in  P  ,  (e.g.,  any 
movement  requirement  can  be  sent  via  any  channel)  then  it  is  appealing  to 
represent  the  aggregate  node  capacities  as  the  sum  of  the  individual  node 
capacities.  This  Implies  that  t^  ■  1. 

The  motivation  for  a  choice  of  t  less  than  1,  is  to  prevent  the 
aggregate  capacity  from  over  estimating  the  actual  capacity  of  the 
system.  For  example,  if  some  movement  requirements  can  only  move  as 
outsize  cargo,  an  aggregate  node  representing  both  oversize  and  outsize 
cargo  would  overestimate  capability  if  it  were  simply  the  sum  of  the 
individual  detailed  oversize  and  outsize  channel  capabilities.  The 
amount  of  usable  capability  can  be  controlled  by  the  factor  t ^ . 

When  node  aggregation  is  performed,  information  is  lost  regarding 
the  detailed  capacities  on  nodes  to  which  each  of  the  multiple  arcs  are 
incident.  The  problem,  as  shown  in  Figure  2.5  represents  the  aggregate 
model  after  the  parallel  arcs  in  Figure  2.4  have  been  aggregated. 

If  a  relationship  between  the  detailed  arc  flows  x^  and  the 
aggregate  arc  flow  y£  is  assumed,  then  disaggregation  is  greatly 
simplified.  A  simple  relation  is 


Cjj  -  s1J  *  yi  where  jePr 
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where  s^  is  the  fraction  of  the  aggregate  arc  flow  which  is  to  be 
allocated  to  the  detailed  arc  lj. 

Such  a  relationship  requires  that  the  values  of  8^  be  specified. 

For  a  given  t^,  if  s^  is  assumed  to  be  a  function  of  the  data  regarding 
sink  nodes  alone,  (i.e.  Independent  of  the  source  node  data)  then  the 
condition  that  x^  be  a  feasible  solution  to  the  detailed  problem  permits 
the  determination  of  the  appropriate  value  for  s^. 

From  the  detailed  problem,  the  aggregate  problem  is  formulated 
corresponding  to  the  network  shown  in  Figure  2.5.  The  aggregate  node 
capacities  are  as  defined  earlier  and  y^  refer  to  the  flows  on  the 
aggregate  arcs  (i,r). 

Aggregate  Problem. 
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If  t  ■  1  is  selected,  and 


8ij  yi 


then,  in  order  to  ensure  that  the  obtained  from  the  solution  to  the 
aggregate  problem  is  feasible  to  the  original  problem,  the  feasibility 
conditions  must  be  checked  as  follows 
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From  the  feasibility  conditions 
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Next,  the  conservation  of  flow  conditions  are  checked  at  the  sink  node 
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Assuming  that  s 
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depends  on  the  sink  data  alone  then 


The  above  relation  holds  for 
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If  8  is  independent  of  the  source  nodes  and  t  «  1  then 
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Also,  for  this  specification  of  t^  and  s^  values,  we  have  the  objective 
function  value  of  the  feasible  solution  generated  to  the  original  problem 
equal  to  the  objective  function  value  of  the  aggregate  problem,  if  the 
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aggregate  arc  costs  are 
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This  follows  since  the  objective  function  of  the  detailed  problem  is 
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for  the  cost  of  the  aggregate  arc  as  specified  above. 


vs  s  S' 


If  the  values  of  the  node  capacities  Involved  in  an  aggregate  node 
are  equal  and  tj  -  1,  then  s^  -  l/(  P^)  for  j  in  Pf.  This  essentially 
means  that  the  arcs  flows  are  distributed  equally  across  all  the  incident 
detailed  arcs. 

The  discussion  in  this  section  permits  the  use  of  a  detailed  problem 
(assuming  that  there  is  an  arc  from  every  source  to  every  single  node) 
together  with  a  specified  aggregation  of  nodes  in  order  to  define  an 
aggregate  problem  with  aggregate  costs 


where 

SiJ 

and  aggregate  node  capacities 


Once  the  aggregate  problem  Is  solved,  a  disaggregation  scheme  is 
available  which  sets 

Xij  *  8ij7i* 

This  always  provides  a  feasible  solution  to  the  original  problem. 
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Unfortunately,  In  many  deployment  planning  situations,  the 
transportation  problem  is  sparse  with  the  sparsity  reflecting  feasible 
allocations  of  movements  to  transportation  modes  (i.e.  some  movement 
requirements  cannot  be  assigned  to  certain  channels).  In  such  cases 
there  is  not  an  arc  from  every  source  node  to  every  sink  node  and  the 
disaggregation  outlined  in  this  section  will  not  always  provide  a 
feasible  solution. 

The  solution  to  the  detailed  problem  generated  is  guaranteed  to  be 
feasible  only  if  all  such  nodes  in  a  cluster  have  arcs  to  the  same  set  of 
source  nodes.  Since  a  specified  aggregation  may  not  satisfy  this 
condition,  formulations  need  to  be  examined  which  generate  feasible 
solutions  when  sink  nodes  with  arcs  to  different  source  nodes  are 
aggregated.  Such  a  procedure  will  be  provided  in  Section  5.0. 
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3.0  DISAGGREGATION  ISSUES: 

Given  a  detailed  problem,  an  aggregate  model  can  be  constructed  as 
Illustrated  In  the  previous  section.  Once  the  aggregate  model  Is  solved, 
the  next  question  Is  hov  to  disaggregate  the  aggregate  solution. 
Aggregation  Is  used  to  determine  a  gross  allocation  of  resources  to  the 
aggregate  sink  nodes.  Disaggregation  distributes  the  allocated  aggregate 
resources  across  the  detailed  nodes  and  arcs. 

Aggregation  is  usually  motivated  by  a  desire  to  decompose  a  large 
problem  into  smaller  manageable  subproblems.  A  desirable  feature  of  any 
disaggregation  is  that  it  permit  the  subproblems  to  be  solved 
independently.  This  permits  the  subproblem  solution  process  to  occur  in 
parallel.  The  objective  of  disaggregation  is  usually  to  generate  "good 
feasible"  solutions  to  the  detailed  problem.  It  is  also  desirable  to  be 
able  to  evaluate  alternative  solutions  to  the  detailed  problems.  The 
characteristics  of  a  good  disaggregation  procedure  are  that  it: 

(1)  enable  parallel  processing  of  the  clusters,  so  that  the  time  is 
merely  the  largest  time  among  all  clusters.  This  is  useful 
from  the  perspective  of  use  of  aggregation  as  a  co-ordinating 
procedure. 

(2)  generate  a  feasible  solution  to  the  original  problem.  This 
permits  the  procedure  to  be  interrupted  by  the  user  interface 
so  as  to  provide  input  to  speed  up  the  overall  problem 
solution.  Also,  even  though  optimality  of  the  entire  problem, 
may  not  have  been  reached  the  solution  at  the  end  of  each 
iteration  offers  suboptimal  scenarios  which  may  be  satisfactory 
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from  the  perspective  of  satisfying  certain  other  qualitative 
considerations . 

Disaggregation  procedures  reported  in  the  literature  are  essentially 
of  two  types:  fixed  weight  disaggregation  and  optimal  disaggregation. 
Fixed  weight  disaggregation  uses  the  same  weights  used  in  the  aggregated 
problem  formulation.  It  is  fixed  in  the  sense  that  it  is  a  single, 
noniterative  multiplication  procedure.  It  uses  the  s  values  to 
generate  a  feasible  solution  to  the  detailed  problem  as  follows. 


‘lj  "  *lj  *  yl  for  JePr 


This  provides  a  feasible  solution  to  the  detailed  problem  provided 
that  there  is  an  arc  from  every  source  node  to  every  sink  node.  This 
procedure  is  easy,  quick,  and  does  not  involve  any  specific  algorithmic 
procedures.  However,  since  the  s^  values  were  assumed  to  be  independent 
of  the  source  nodes,  the  solution  generated  for  sparse  problems  may  be 
infeasible  because  information  regarding  the  missing  arcs  is  not 
utilized . 

Optimal  disaggregation  is  optimal  in  the  sense  that,  given  a 
partition  of  the  source  flow  allocation  to  the  clusters,  this  flow  Is 
routed  optimally  within  each  cluster.  The  solution  generated  is  no  worse 
then  the  fixed  disaggregation  case. 

The  optimal  disaggregation  procedure  allows  each  cluster  of  nodes  to 
be  evaluated  independently  from  the  aggregate  problem  solution.  This 
permits  parallel  solution  of  the  disaggregation  problems  for  each 
cluster.  The  essential  idea  is  to  flow  the  aggregate  flows  back 
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through  the  detailed  arcs  and  nodes.  The  problem  constructed  for  each 
cluster  would  be  a  transportation  problem.  It  would  have  a  number  of 
source  nodes  equal  to  the  number  of  aggregate  basic  feasible  arcs 
incident  to  the  aggregate  cluster  node,  and  a  number  of  sink  nodes  equal 
to  the  number  of  detailed  nodes  aggregated  into  the  cluster.  The  costs 
on  the  arcs  are  the  costs  on  the  detailed  problem,  with  the  aggregate 
flows  forming  the  capacities  on  the  source  nodes  and  capacities  (obtained 
from  the  detailed  problem)  on  the  sink  nodes.  The  disaggregation  problem 
for  each  cluster  r  is  as  follows 
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jePr  J 

for 

i  -  1,2, ...M 

S  x  *  b 

l  ij  j 

for 

J£Pr 

X 

p* 

L-t. 

V 

o 

The  corresponding  general  network  model  is  shown  in  Figure  3.1. 

Since  the  disaggregation  problems  allow  flows  only  on  arcs  which  exist  in 
the  original  problem,  a  feasible  solution,  if  one  exists  is  assured  to 
the  detailed  problem. 

For  the  example  whose  detailed  formulation  is  given  in  (I)  and  whose 
aggregate  formulation  is  given  in  (II),  disaggregation  involves 
allocating  the  aggregate  flow  variables  y^  to  the  detailed  flow  variables 
x^j.  Since  there  are  two  aggregate  clusters,  there  are  two  disagregation 
problems . 


19 


For  each  cluster,  the  following  Is  the  optimal  disaggregation 
problem.  The  values  of  are  specified  from  the  solution  to  the 
aggregate  problem.  The  networks  are  shown  in  Figure  3.2  for  cluster  1 
and  in  Figure  3.3  for  cluster  2.  The  constraints  corresponding  to  these 
networks  are  given  in  (III)  and  (IV)  respectively. 


(13)  xu  +  x15 


(16)  x. 


+  X25 


dS)  ^ 


+  x36  +  x3?  -  y3 


(III) 


(21)  x. 


Si 


4.0  THE  ITERATIVE  STEP 

Once  an  aggregate  problem  has  been  solved  and  a  feasible  solution  to 
the  detailed  problem  generated,  this  Information  Is  used  to  decide  how 
the  aggregate  problem  can  to  be  modified  so  that  the  procedure  will  move 
towards  an  optimal  solution  of  the  detailed  problem.  The  goal  is  to  add 
enough  additional  information  to  the  aggregate  problem  to  move  It 
ultimately  to  an  optimal  solution  of  the  detailed  problem. 

Methods  have  been  developed  to  generate  a  dual  variable  solution  to 
the  detailed  problem  using  the  aggregate  dual  information.  These  dual 
variable  values  can  be  used  to  identify  dual  infeasible  arcs  in  the 
detailed  problem.  (Dual  infeasible  arcs  are  those  which  do  not  satisfy 
the  dual  constraints.)  Given  dual  infeasible  arcs,  various  procedures 
have  been  suggested  to  reaggregate  the  detailed  problem.  Balas  [1] 
suggests  the  following  procedure.  The  aggregate  problem  duals  are  used 
to  generate  a  dual  variable  solution  to  the  detailed  problem.  All  the 
nodes  are  separated  into  a  cluster  that  have  dual  infeasible  arcs 
incident  to  them.  Singleton  clusters  are  setup  for  each  such  node  in  the 
new  aggregation.  This  procedure  is  iterated  until  the  solution  of  the 
detailed  problem  is  obtained. 

The  problem  with  this  procedure  is  that  the  aggregate  problem  size 
tends  to  grow  very  rapidly.  This  may  not  be  a  desirable  feature. 

Taylor  [4]  suggests  a  modified  approach  where  only  the  most  dual 
infeasible  arc  in  each  cluster  is  considered.  The  corresponding  node  tc 
which  it  is  incident  is  separated  to  form  a  single  node  cluster.  This 
method  uses  the  same  motivation  as  the  earlier  approach  but  controls  the 
growth  of  the  aggregate  problem  size.  The  procedure  uses  the  duals  to 
identify  the  most  dual  infeasible  arc  in  a  cluster.  From  each  cluster  r. 
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Che  node  j^eP^  is  separated  Co  form  a  cluster  by  Itself  in  the 
reaggregated  problem. 

These  procedures  can  be  summarized  as  follows: 

(1)  The  problem  Information  is  used  to  partition  nodes  into 
clusters  of  similar  nodes. 

(2)  These  nodes  are  aggregated  to  setup  an  aggregate  problem 

(3)  The  aggregate  problem  solution  is  used  to  generate  upper  and 
lower  bounds  on  the  original  problem. 

(4)  The  dual  aggregate  solution  is  used  to  generate  a  dual 
solution.  This  dual  solution  is  used  along  with  additional 
information  about  the  problem  instance,  to  reaggregate  the 
problem. 

(5)  Steps  (1)  through  (4)  are  repeated  until  an  optimal  solution  to 
the  detailed  problem  is  obtained  or  until  the  user  is  satisfied 
with  the  sub-optimal  solutions  generated  at  some  stage. 

These  procedures  have  two  major  shortcomings  with  regard  to 
deployment  planning  problems.  The  first  is  that  they  only  work  for  the 
situation  where  energy  node  in  a  cluster  has  the  same  connectivity  to 
source  nodes  (i.e.,  if  one  sink  and  in  the  cluster  is  connected  to  a 
particular  source  node,  then  all  sink  nodes  in  '.he  cluster  are  connected 
to  the  same  source  node).  This  condition  is  not  normally  satisfied  for 
the  MR  assignment  problem. 

The  second  shortcoming  is  that  a  larger  aggregate  problem  results 
after  each  iteration.  If  the  computer  system  is  constrained  by  storage 
space,  then  this  can  result  ultimately  in  an  unacceptably  large  problem. 

The  following  reformulation  has  the  potential  to  over  come  both  of 
these  problems. 
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5.0  MODEL  REFORMULATION 


For  the  example  problem  in  section  2.0,  consider  the  original 
constraints  (1)  through  (7)  and  the  new  constraints  (8)  through  (22). 
Substituting  constraints  (13), (14),  and  (15),  constraint  (11)  can  be 
rewritten  as 


X14  +  X15  +  X25  +  X34  “  a4  +  a5 


This  is  a  relaxation  of  the  original  constraints  (4)  and  (5).  Similarly 
constraint  (12)  is  a  relaxation  of  constraints  (6)  and  (7). 

Constraint  (8)  together  with  constraints  (13)  and  (18)  provides 


X14  +  X15  +  X16  "  al 


which  is  constraint  (1).  Constraint  (9)  with  (14),  and  (19)  provides 
(2);  constraints  (10)  with  (15)  and  (20)  provides  (3).  Constraints  (4) 
through  (7)  are  Identical  to  (16)  through  (17)  and  (21)  through  (22). 
Thus,  if  constraints  (8)  -  (22)  were  compared  to  the  original  problem  (1) 
through  (7),  the  additional  constraints  are  simply  relaxations  or  can  be 
used  to  generate  the  same  original  constraints.  Constraint  set  (8) 
through  (22)  can  be  rearranged  to  form  the  system  of  equations 


Min  I  I  c,  .  x 

1  j 

subject  to:  constraints  (8)  through  (12) 
and  for  cluster  1  :  constraints  (13)  through  (17) 
for  cluster  2  :  constraints  (18)  through  (22) 


This  system  of  equations  has  a  matrix  structure  as  shown  in  Figure 


5.1.  The  matrix  has  a  structure  which  consists  of  a  master 
transportation  problem  linked  via  columns  to  subproblems  which  are  also 
transportation  problems. 


5.1  General  Model 

The  above  specific  example  can  be  generalized  with  the  following 
notation 


Detailed  Problem 

:  Min  lie  x 

ij  J  3 

(23) 

s.t.  Z  x  -  a 

j  3 

i  -  1...M 

(VI) 

(24) 

l  ‘  bj 

j  -  1...N 

'ij  >0 

Given  a  specified  aggregation  of  detailed  nodes  an  equivalent 


problem  can  be  formulated  as  follows. 

Let  r  *  1,2,... L  be  the  number  of  clusters. 

Let  denote  the  nodes  j  which  comprise  each  cluster  r. 

Let  y^  be  defined  as 

(25)  y^  -  E  x.,  for  each  i  ■  1,2, ...M  and  r  *  1,2, ...L  For 

1  eP 
J  r 

each  r  -  1,2,...L  add  the  redundant  constraints 
M 

(26)  Z  y!  -  Z  b 

1-1  j£Pr  3 

to  problem  (V).  Constraints  (26)  are  redundant  since  combined  with  the 
definitions  of  y^  In  (25)  a  relaxation  of  the  detailed  problem 
constraints  (24)  is  obtained.  The  y^  variable  definitions  (25)  also 
result  in  an  equivalent  representation  of  the  constraints  (23)  as 
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27 


L 

(27)  £  y!  -  a±  1-1,2,..* 

T-\ 

The  constraint  sets  (24)  through  (27)  can  be  rearranged  to  obtain 


Hin 


£  £  c 


ij  U 


subject  to:  £  y*  -  a  i  *  1,2,.. M 

i-1 


M 

£  yT.  -  £  b  r  -  1,2,  ..L 

i-1  je?r  3 

(VI) 

and  for  each  cluster  r,  r  -  1,2,..L 

E  Xi1  “  yi  1  " 

j£Pr  J 


l  "  bJ  jEPr 


x 


ij’ 


y\  >  o 


The  matrix  structure  for  the  detailed  problem  is  shown  in  Figure  5.2,  and 
that  for  the  reformulation  is  shown  in  Figure  5.3. 

Every  feasible  set  of  x^  for  (VI)  is  also  feasible  to  (V),  and  vice 
versa.  This  follows  since  only  redundent  constraints  have  been  added  to 
(V)  in  order  to  obtain  (VI).  Hence,  solving  (VI)  also  solves  (V).  A 
reformulation  of  the  form  (VI)  can  be  constructed  for  any  specified 
aggregation  of  nodes. 
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V 


5.2  Solution  of  the  Reformulated  problem 

There  are  two  extreme  cases  with  regard  to  the  reformulated 
problem.  If  each  node  forms  a  cluster  then  the  aggregate  problem  Is  as 
large  as  the  original  problem.  If  all  nodes  together  form  a  cluster 
then  the  aggregate  problem  Is  trivial.  In  the  latter  case  there  Is  just 
one  arc  out  of  each  source  node  to  the  aggregate  sink  node,  and  the 
solution  Is  to  fill  the  arc  to  the  source  supply  value.  Between  these 
two  extremes,  there  Is  a  distribution  of  time  and  complexity  of  solution. 
These  vary  between  the  aggregate  and  detailed  levels  which  allows  the 
problem  to  be  solved  within  the  time  and  space  specifications. 
Considerable  latitude  Is  available  In  tradlng-off  size  the  aggregate 
problem  versus  size  and  number  of  the  disaggregated  problems. 

An  immediate  approach  to  solving  the  reformulated  problem  Is  to 
setup  the  aggregate  problem  as  a  master  problem.  The  subproblems  are 
formed  by  the  disaggregation  problems  linked  together  by  Benders' 
constraints  [2].  For  the  problem  discussed  In  section  5.0,  the  master 
problem  would  be 


subject  to:  ►  yT  -  a, 
r-1  1 


1-1... 3 


E  yi  “  a4  +  a5 
l-l  1 


yi  “  a6  +  a7 


one  Benders'  constraint  for  each  Iteration  and 


y^  >  0  for  1-1. ..2,  r-1. ..2 


rat 
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The  subproblems  for  a  given  y^  would  be 
SUBPROBLEM  1 


Min  cu*u  +  c15x15  +  c25x25  +  c34x34 


*14  +  *15 


*25 


14 


34 


34 


X15  +  *25 


SUBPROBLEM  2 

Min  c16x16  +  <^7X27  +  <^5*36  +  c37*37 
*16 

*27 


*36  +  *37 


16 


+  x 


27 


2 

y2 


36  6 

+  *37  "  “7 

and  x^j  >  0 


Benders'  constraints  of  the  fora 


^  _  1  1_2  2  1  .1  2  2 
z  >  Z  Uf  7i  +  lUi  Xf  +  Vla4  +  v2a5  +  v^  +  v2a? 


are  added  to  the  master  problem  at  each  iteration,  with  the  u’s  and  v* s 
being  the  duals  of  subproblems  1  and  2. 

Benders'  procedure  is  discussed  in  detail  in  PDRC  reports  84-09  and 


85-03.  The  procedure  involves  solving  the  master  problem  (in  this  case 
the  aggregate  problem)  and  then  solving  the  subproblems  (in  this  case  the 


disaggregation  problem).  Based  on  the  solutions  to  the  subproblems,  a 
Benders'  constraint  is  constructed  and  added  to  the  master  problem.  The 
new  master  problem  is  solved  and  the  process  repeated.  A  lower  and  upper 
bound  on  the  value  of  the  optimal  objective  function  of  original  detailed 
problem  is  generated  at  each  iteration.  When  the  two  bounds  are  equal, 
the  process  terminates  with  an  optimal  solution. 


5.3  Advantages  and  Disadvantages  of  the  Benders'  Approach 

Benders'  decomposition  is  suggested  to  solve  the  reformulated 
problems.  The  basic  advantage  of  Benders'  procedure  is  that  it  separates 
the  solution  of  the  aggregate  problem  and  the  disaggregate  subproblems. 
Once  the  aggregate  problem  is  solved,  the  disaggregation  subproblems,  one 
for  each  cluster,  can  be  solved  Independently.  The  linking  Benders' 
constraints  generated  by  the  subproblem,  one  constraint  for  each  cluster, 
are  added  to  the  aggregate  problem,  guiding  it  towards  the  solution  of 
the  entire  problem.  Furthermore,  at  each  iteration  upper  and  lower 
bounds  on  the  optimal  objective  function  value  are  generated.  These 
bounds  permit  termination  when  the  solution  is  reasonably  close  to 
optimal. 

A  problem  with  the  use  of  Benders'  decomposition  is  that  the 
constraints  added  to  the  aggregate  problem  cause  it  to  be  a  network  flow 
with  side  constraints  model.  This  reduces  the  advantage  of  its  original 
pure  network  structure.  Research  is  in  progress  on  a  procedure  which 
permits  the  additional  Benders'  constraints  to  be  used  to  change  flow 
capacities  in  the  aggregate  transportation  model.  This  would  allow  use 
of  network  algorithms  to  solve  the  aggregate  problem.  Based  on 


information  available  at  an  iteration,  the  aggregation  itself  may  also  be 
changed  to  hasten  the  solution. 

While  Benders'  decomposition  is  a  procedure  which  can  be  used  to 
solve  the  reformulated  model,  other  procedures  which  exploit  the  network 
structure  of  the  aggregate  and  disaggregation  problems  many  expedite  the 
solution  mechanism  for  the  reformulated  model. 


6.0  EXTENSIONS  AND  RESEARCH  AREAS 

The  concepts  discussed  In  this  report  suggest  various  areas  for  future 
research. 

(1)  Though  the  aggregation  has  been  discussed  with  respect  to  a  two 
level  solution  procedure,  the  procedure  could  be  extended  to  a 
multi-level  hierarchical  procedure.  The  decisions  are  refined 
at  each  level  with  preceding  level  providing  resource 
constraints  for  succeeding  levels. 

(2)  Since  any  minimus  cost  flow  problem  on  a  general  network  can  be 
transformed  Into  a  transportation  problem  as  In  Lawler  [3],  the 
aggregation  procedure  could  be  applied  to  solve  minimal  cost 
flow  problems.  This  aggregation  and  problem  reformulation 
might  be  carried  out  without  transforming  the  network,  thereby 
working  directly  with  the  network  itself. 

(3)  Since  the  master  and  subproblems  are  transportation  problems, 
techniques  which  replace  the  Bender's  cuts  with  equivalent 
capacities  on  aggregate  arcs  might  be  devised  so  that  the 
master  problem  remains  a  network  flow  problem,  rather  than 
become  a  network  flow  with  side  constraints  model. 

(4)  In  solving  the  mas'.er  problem,  the  Benders'  might  retained  in 
the  objective  function.  This  would  produce  a  network  flow 
problem  with  a  mlnimax  objective  function.  In  this  case  the 
objective  function  is  a  piecewise  linear  convex  objective 
function. 

(5)  Even  though  only  sink  node  aggregation  has  been  discussed,  it 
can  be  extended  to  aggregation  of  both  source  and  sink  nodes  by 
performing  the  operation  sequentially. 
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